Phase Separation Based on U(l) Slave-boson 
Functional Integral Approach to the t-J Model 
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We investigate the phase diagram of phase separation for the hole-doped two dimensional sys- 
tem of antiferromagnetically correlated electrons based on the U(l) slave-boson functional integral 
approach to the t-J model. We show that the phase separation occurs for all values of J/t, that 
is, whether < J/t < 1 or J/t > 1 with J, the Heisenberg coupling constant and t, the hopping 
strength. This is consistent with other numerical studies of hole-doped two dimensional antiferro- 
magnets. The phase separation in the physically interesting J region, < J/t < 0.4 is examined by 
introducing hole-hole (holon-holon) repulsive interaction. We find from this study that with high 
repulsive interaction between holes the phase separation boundary tends to remain robust in this low 
J region, while in the high J region, J/t > 0.4, the phase separation boundary tends to disappear. 



One of the most interesting observations in high- 
Tc cuprates (superconductors) is the phase separation, 
which may play an important role on superconductivity. 
The phase separation results from a thermodynamic in- 
stability which arises from the violation of the stability 
condition, = r?d'^e/dr? — r?d[i/dn > 0. Here K is 
the compressibility; e, the ground state energy per site; 
n, the electron density; and ^, the chemical potential. 
Initially the phase separation instability was believed to 
inhibit superconductivity. Recently it draws a great at- 
tention owing to its possible connection with supercon- 
ductivity based on experimental observations 
in high- Tc cuprate oxides. 

We write the t-J Hamiltonian for the study of the hole 
doped systems of antiferromagnetically correlated elec- 
trons. 
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with = l/2c^^(Ta/3Cj/3 and Ui = J2cr4a'^i<^ where c}^ 
creates an electron of spin a on site i. Earlier, using the 
t-J model a possibility of phase separation in high-Tc 
cuprates has been brought up by Emery et al. [||. They 
predicted the existence of phase separation at all possible 
values of J/t, that is, < J/t < 1 or J/t > 1 where J is 
the antiferromagnetic correlation strength and t, the hop- 
ping integral, including the case of J/t < 1. On the other 
hand, other numerical studies ||(HfOfl predicted the exis- 
tence of phase separation only for J/t > 1 where J value 
is unrealistic for the high Tc cuprates of current inter- 
est. Recently, from a Green function Monte Carlo study 
Hellberg and Manousakis reported that the phase 
separation can occur for all values of J, in agreement 
with the earlier exact diagonalization study of Emery et 
al Q. In the present study, by using the U(l) slave-boson 
functional integral method |12-17|, we obtain a phase di- 
agram in the plane of electron density vs. J/t, by using 
the Maxwell construction pjisl. 



If violation of the stability condition K^^ > occurs 
in the electron density range of rti < Tie < n2, where 
rii is the electron density for a hole- rich phase and n2, 
the electron density for a hole-free phase, the system is 
expected to separate into two subsystems with electron 
densities ni and n2 respectively. Since we are interested 
in the hole-doped systems of high Tc cuprate oxides, the 
physics can be conveniently described in terms of the 
hole density, x = 1 — rig. Thus we first examine the 
ground state energy density eh{x) of the hole doped sys- 
tem as a function of hole density x. For the Maxwell's 
construction [p|p^ to treat a finite system, we consider a 
straight line which intercepts both the Heisenberg energy 
en — Sh{0) at x = and a curve given by eh{x) ai x ^ 0. 
The slope of the straight line at x is then given by 



e{x) 
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If a minimum of e{x) is found at a hole density of x = Xc 
(at which the straight line intercepts tangentially the 
curve of eh{x)), phase separation is expected to exist be- 
low the onset (critical) density of x = a;c (sjjl^. As a 
result the energy density of the phase separated system 
is described by the the linear function with slope of e{xc) 
in the doping range of < a; < Xc. Thus in this region 
the system is stabilized with its energy lower than that 
of the uniform phase, by forming a system composed of 
two subsystems: one with a hole-rich phase of the elec- 
tron density of ni = 1 — Xc and the other with a hole-free 
phase of the electron density of n2 = 1. 

In order to clarify how to compute the ground state 
energy as a function of electron or hole density we briefly 
discuss our earlier approach |17| of the U(l) slave-boson 
representation of the t-J Hamiltonian. In this approach 
we introduce an additional contribution of hole-hole re- 
pulsion to the original t-J Hamiltonian, 
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with = l/2f}^aai3fi(3 and rn = Y.a fLfi<y A^o is 
the chemical potential to fix the number of electron 
to Ne- Here the local constraint of single occupancy, 
J2a fLf^T + = 1 is assumed, flifia) is the spinon 
creation (annihilation) operator and the holon 

annihilation (creation) operator. The nearest neighbor 
(NN) configuration of two holes is energetically more fa- 
vorable than other possible configurations. This is evi- 
dent from the separate inspection of the two attractive 
interaction terms, JS^-Sj and —(J/4) riirij. For the latter 
we write IITtI , 
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where the local constraint of single occupancy is taken 



into account. The effective attraction between the NN 
holes arises from the last term of the equation above. In 
view of the numerical finding of excessive large binding 
of hole pairs |^ , the hole-hole repulsive interaction term 
(the third term in Eq. (||)) is introduced. Owing to its 
introduction, we are now able to examine how phase di- 
agram boundary is affected by the variation of hole-hole 
repulsion interaction V. Such discussion will be made 
later. 

As a result of the Hubbard-Stratonovich transforma- 
tion [|l6| in Eq. (||), the Heisenberg exchange and the 
hopping terms are led to linearized terms involving the 
hopping order field Xji = {St/SJb^^bi + f]„fia) in asso- 
ciation with the exchange interaction channel and the 
spinon singlet pairing order field AJ^ = (/jt/a - /ji/it) 
in association with the pairing channel. The contribu- 
tion of the direct (Hatree) channel is omitted based on 
the assumption of paramagnetic states for each site, i.e., 
(Si) =0 ijl^. Long-range antiferromagnetic fluctuations 
are thus ignored in this approach Jl6| , pTt . The resulting 
effective Hamiltonian is then , 
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where a Lagrange multiplier field is introduced to im- 
pose the local constraint of single occupancy for both the 
spinon and the holon. will be absorbed into the effec- 
tive chemical potential of holon and spinon in the mean 
field approaches . The quartic holon term (the second 



term in Eq. (||) above) ^ j) i^ repulsive 
pO| . It is important to realize that this term involves 
nothing but the holon exchange interaction. Thus al- 
lowing the holon exchange channel for the quartic holon 
term (the second term in Eq. (|^)), we obtain 
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We find from our numerical calculation of the Maxwell construction that the above holon exchange channel also 
affects phase separation, by effectively reducing the hopping (kinetic) energy. We find that with the neglect of the 
contribution of Eq. (^) the phase separation does not occur even at sufficiently low doping. 

The effective holon attractive interaction term (when < < J/4 in the third term of Eq. (||)) ~ 

V)blb'jbibj can be decomposed into terms involving the direct, exchange, and pairing channels. The Hubbard- 
Stratonovich transformations and Bogoliubov-Valantin transformation in the momentum-space are made. As de- 
scribed in Ref. Il^], the mean field free energy at hole doping rate x is then. 
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Here E: 
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" l^-^y + is the quasi-particle exci 

tation energy for spinons and E'^ = \J {e\ — fi'^Y ~ 

for holons, where and /i'' are the effective chem- 
ical potentials of spinon and holon respectively and 
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^Xlk, 4 = -2icffX7fc with 7^ = cosfca; 



cos ky . The effective Heisenberg coupling constant is 
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ping strength is teff = 2^*+8V(3j) 
of spinon and holon order parameter 77 



and the effective holon hop- 

V 



where r/ is the ratio 
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From the minimization of the mean field free energy 
with respect to the scalar fields, x, A^, and A**, we deter- 
mine the ground state energy density of the hole-doped 
system per site, eh{x) = \m\T^o Fmf{Xi ^ ^'') / ^ as a 
function of hole density x. The critical hole density is 
actually found from the Maxwell construction by intro- 
ducing e(x) = {eh{x) — e}i)lx. In Fig. 1 we display the 
Maxwell's construction for J = l.ht and F = 0. The in- 
set is the Green function Monte Carlos calculation of the 
t-J model by Hellberg and Manousakis jl^]. Quantitative 
disagreement exists between our U(l) slave-boson func- 
tional integral approach and the Green function Monte 
Carlo calculation. However we find that there exists a 
minimum of e(x) at the critical hole density of x^ 0.72 
which is close to the value of Hellberg and Manousakis 
Iprf , Xc ~ 0.7, below which phase separation occurs. 

In Fig. 2 we display a predicted phase diagram in 
the plane of the unitless Heisenberg exchange coupling 
strength, J jt and the electron density, n = \ — x for 
various values of NN hole-hole (holon-holon) repulsion 
energy V . The phase separation was predicted to occur 
for J/t < 1 and J/t > 1. This prediction is consis- 
tent with other numerical studies Green function 
Monte Carlo resuhs (sohd hne) up to the 28 x 28 
square lattice and and the exact diagonalization result 
[^1 (stars) with a 4 x 4 lattice a are displayed for compar- 
ison with our results (solid circles for V = 0) with 




hole density x = 1 - 

FIG. 1. Maxwell's construction: e{x) vs. x for J = O.lt, 
V — 0. The Green function Monte Carlo result [11] in the 
inset is denoted by HM. 




FIG. 2. Phase separation for the hole-doped systems 
of antiferromagnetically correlated electron in the plane of 
Heisenberg coupling strength, J/t and the electron density, 
n = 1 — X. The solid line denoted by HM is the Monte Carlo 
prediction (Ref. [11]); and the stars denoted by EKL, the 
result of Emery et al (Ref. [5]); and the solid circles are our 
computed results for V = 0. The phase separation boundaries 
for V = 0.125 J (open circles) and V = 0.2 J (open triangles) 
are also displayed. The critical hole doping density Xc is seen 
to decrease with the increase of V in the region of large J/t. 



100 X 100 lattice [^. Despite some numerical differ- 
ences, interestingly all of these methods yield nearly the 
same critical Jc ~ 3.4t, above which the hole-rich phase 
contains no spins, as shown in the Fig. 2. In the small 
J/t limit, the phase separation is expected to occur as 
a result of relative increase in kinetic energy (compared 
to the Heisenberg interaction energy J) which promotes 
relatively easier hopping of holon (holes) from site to 
site, thus avoiding antiferromagnetic spinon (spin) frus- 
trations to lower the energy of the system and creating a 
hole-rich region. On the other hand, in the large J limit 
the phase separation occurs owing to the Heisenberg in- 
teraction coupling which promotes the antiferromagnetic 
order by inhibiting the occurrence of holes in the region 
of the antiferromagnetic phase 

We now explore the effects of hole-hole (holon-holon) 
repulsive interaction V on the phase separation bound- 
ary. The uniform phase is expected to be more favorable 
owing to the enhanced difficulty of hole pairing with the 
increase of V. The critical doping density in the phys- 
ically interesting J region, J/t < 0.4 is predicted to be 
relatively insensitive to the variation of V compared to 
the high J limit, as shown in Fig. 2. We observe that the 
critical hole doping density Xc quickly decrease beyond 
the large J region of J/t > 0.4 as the holon-holon (hole- 
hole) repulsion energy V increases. For the case of the 
large V limit {V ~ 0.25 J), the uniform phase occurs with 
a small critical hole density, indicating the phase separa- 
tion is not likely to occur. As shown in Fig. 2, we note 
the persistence of phase separation in the region of small 
J and the propensity of gradual disappearance of phase 
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separation in the region of large J /t at this high Umit of 
the hole-hole repulsive interaction. Such persistence of 
phase separation at small J/ 1 despite the increase of V is 
attributed to the effective increase of the kinetic energy 
of holons, to avoid the frustration of antiferromagnetic 
spinous g]. 

In the present study we investigated the phase diagram 
involving phase separation based on the U(l) slave-boson 
functional integral approach to the t-J model. We find 
that the phase separation occurs in the region of low hole 
doping for all possible values of Heisenberg coupling con- 
stant J, that is, whether J /t < 1 or J /t > 1 with an 
upper bound of J/t (at J ~ 4.2t). This observation is 
consistent with other numerical studies of hole-doped two 
dimensional antifcrromagnets. 



[1] C. Castellani, C. Di Castro, and M. Grilli, 

Phys. Rev. Lett. 75, 4650 (1995). 
[2] V. J. Emery, S. A. Kivelson, and O. Zachar, Phys. Rev. B 

56, 6120 (1997). 
[3] J. M. Tranquada, J. Phys. Chem. Solids 59, 2150 (1998); 

references therein. 
[4] A. W. Hunt, P. M. Singer, K. R. Thurber, and T. Imai, 

Phys. Rev. Lett. 82, 4300 (1999); references therein. 
[5] V. J. Emery, S. A. Kivelson, and H. Q. Lin, 



Phys. Rev. Lett. 64, 475 (1990). 
[6] H. Fehske, V. Waas, V. Roder, and H. Biittner, 

Phys. Rev. B 44, 8473 (1991). 
[7] W. O. Putikka, M. U. Luchini, and T. M. Rice, 

Phys. Rev. Lett. 68, 538 (1992). 
[8] D. Poilblanc, Phys. Rev. B 52, 9201 (1995). 
[9] H. Yokoyama and M. Ogata, J. Phys. Soc. Janpan, 65, 

3615 (1996). 
[10] M. Kohno, Phys. Rev. B 55, 1435 (1997). 
[11] C. S. Hellberg and E. Manousakis, Phys. Rev. Lett. 78, 

4609 (1997). 

[12] G. Kotliar and J. Liu, Phys. Rev. B 38, 5142 (1988). 

[13] Y. Suzumura, Y. Hasegawa, and H. Fukuyama, 
J. Phys. Soc. Jpn. 57, 401 (1988); H. Fukuyama, 
Prog. Theo. Phys. Suppl. 108, 287 (1992). 

[14] L. B. loffe and A. L Larkin, Phys. Rev. B 39, 8988 (1989). 

[15] P. A. Lee and N. Nagaosa, Phys. Rev. B 46, 5621 (1992). 

[16] M. U. Ubbens and P. A. Lee, Phys. Rev. B 46, 8434 
(1992). 

[17] T. -H. Gimm, S. -S. Lee, S. -P. Hong, and S. -H. S. Salk, 

Phys. Rev. B 60, 6324 (1999). 
[18] A. C. Cosentini, M. Capone, L. Guidoni, G. B. Bachelet, 

Phys. Rev. B 58, 14685 (1998); references therein. 
[19] C. Gazza, G. B. Martins, J. Riera, and E. Dagotto, 

Phys. Rev. B 59, 709 (1999). 
[20] G. Baskaran, Phys. Scr. T 27, 53 (1989). 
[21] Our numerical results are based on Eq. (Q) at T = 

limit which is expressed in momentum space. Numerical 

convergence was achieved for a 100 x 100 square lattice, 

which yielded a good convergence. 



4 



